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Abstract. We suggest that the presence of a quantum gravity induced minimal 
length can be explored using neutrino oscillation probabilities. Neutrinos seem ideally 
suited for this investigation because they can propagate freely over large distances 
and can therefore pile up minimal length effects beyond detectable thresholds. We 
determine the modified survival probability in a scenario with a minimal length and 
find deviations from the classical behaviour for high energies. We find that for the 
currently available experimental statistics the deviations from the standard oscillations 
do only allow for a bound of £ _1 > 10 GeV from MINOS data. On the other hand, 
oscillations of high-energy neutrinos emitted by galactic and extragalactic sources are 
strongly suppressed, leading to a possible observation of quantum gravity effects at 
neutrino telescopes such as IceCube and ANTARES. 
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1. Introduction 

Quantum gravity has been drawing the attention of the scientific community for more 
than thirty years, but it is still far from being fully understood and explained. While 
there is some progress on the theoretical side due to the presence of at least two 
viable candidate theories, string theory and loop quantum gravity, the absence of any 
experimental quantum gravity data up to now is still puzzling [1]. The main difficulty 
in getting experimental evidence is connected to the extremely large value of the Planck 
energy scale Ep ~ 10 16 TeV at which quantum gravity signals are supposed to appear. 
This is 15 orders of magnitude higher than the energy scale tested by the LHC. However, 
the arguments leading to this estimation rely on the assumption that general relativity 
holds up to the Planck scale. New physics might enter much earlier and make quantum 
gravity effects important at much lower energy scales. 

In recent years a new field has emerged that does not aim at a fundamental description 
of quantum gravity but instead uses effective theories to study quantum gravity effects 
phenomenologically. Typically, one starts from the conventional theories, general 
relativity or quantum field theory, and implements an effect that is motivated by a 
fundamental theory. A nice overview over the field of quantum gravity phenomenology 
can be found in [2]. 

The effect we will concern ourselves with is the emergence of a minimal length. This 
idea is dated back to early times of quantum gravity [3] and emerges in all formulations 
of quantum gravity (cf. [4], [5], [6]). Different approaches to implement a minimal 
length exist and give rise to a model of quantum spacetime. For instance, generalised 
uncertainty principle (GUP) approaches [7], [8] start from a modified commutator 



and work in a momentum representation of the corresponding algebra. In this paper, 
we will work with a non- commutative geometry (NCG) approach (for a review of other 
approaches see [9]). In [10] an original formulation based on coordinate coherent state 
NCG has been successfully employed to improve classical curvature singularities in black 
hole spacetimes [11, 12] and provide a reliable description of the fractal structure of the 
universe at the fundamental scale [13]. Starting from a commutator 



a representation of the corresponding algebra in terms of coherent position states is built. 
In this representation, delta functions appearing in particle field theories are widened 
into Gaussian distributions. Equivalently to starting from the coherent states, it has 
been shown that the non-commutative smearing effect can be represented by the action 
of a non-local operator e f2Ax which spreads the point-like Dirac delta into a Gaussian 
distribution, as well, i.e. 



[x,p] = t(l + f(p 2 )) 
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where A x is the Laplacian operator onad dimensional Euclidean manifold and £ is the 
NCG induced minimal length. In momentum space, this corresponds to modifying the 
integration measure as 

d 3 p -> e- e2p2 d 3 p. (4) 

The leading order corrections to any field equation in the presence of a non-commutative 
background can then be obtained by replacing the conventional point like source 
term (matter sector) with a Gaussian distribution, while keeping differential operators 
(geometry sector) formally unchanged. The main approach to non-commutative 
geometry via the Moyal ★-product suffers from the necessity to expand scattering 
amplitudes in the minimal length parameter, destroying the non-local character of the 
theory. In contrast, our deformation in Eq.(4) contains an infinite number of terms and 
is intrinsically non-local, leading to a UV finite and unitary field theory [10, 11]. 
In the above formulation the minimal length is not set a priori. Typically, one expects i 
to be of the order of the Planck length. However, as mentioned above, quantum gravity 
effects might become stronger at much lower scales. All we know from experiments 
is that i > 1 TeV [14]. In this paper, we propose to study neutrino oscillations as 
an accurate phenomenon to test quantum gravity effects resulting from the minimal 
length. A modified oscillation pattern could then be used to constrain the value of the 
minimal length. Note again that this is a phenomenological study. We do not provide a 
quantitative analysis but aim for a qualitative study of the minimal length effects and 
their possible observation. 

High-energetic particles probe the structure of spacetime to microscopic scales and 
should be able to feel effects induced by a minimal length. For this reason neutrino 
propagation has already been the subject of investigations for possible tests of quantum 
decoherence [15], modified dispersion relations [16] and modified de Broglie relations [17]. 
In addition, quantum gravity effects in neutrino physics have also been investigated 
by the IceCube collaboration [18]. In contrast to these studies, we study neutrino 
oscillations in a framework that has already been successfully employed to study minimal 
length effects in collider experiments [19], the Casimir effect [20], the Unruh effect [21], 
the magnetic moment of the muon [22] and in black hole physics [11, 12]. 
Even if the expected deviations are very small, neutrinos seem ideally suited for minimal 
length studies. Neutrinos can propagate freely over astronomic distances without 
interacting with matter, giving rise to the hope that the propagation over long distances 
will pile up the minimal length effect beyond currently detectable thresholds. As a simple 
analogy from mechanics, friction effects of a moving particle become significant for high 
velocities, but also for small velocities if the propagation distance is long enough. 

2. Neutrino oscillations 

The most common interpretation of neutrino oscillations is that neutrinos do not 
propagate in a flavour eigenstate but in a mass eigenstate. Neutrinos are not massless, as 
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assumed in the Standard Model, but have a non- vanishing mass. The basis change from 
the flavour eigenbasis to the mass eigenbasis is described by a unitary 3x3 -matrix U 
(the Pontecorvo-Maki-Nakagawa-Sakata matrix) which is parametrized by three mixing 
angles 6 , 12 , 9 i3 , #23 and a CP- violating phase 5q p- 

3 

Wa) = J2 U ^Wk), (5) 
k=l 

where Greek indices stand for flavour eigenstates \u a ) while Latin indices stand for mass 
eigenstates \u k ). Since the free Hamilton operator is diagonal in the mass eigenbasis we 
have 

\u k (t)) = exp(-iE k t) \u k ) , (6) 

with E k being the energy of the neutrino and t the propagation time. The oscillation 
propability for a flavour change from flavour a to flavour j3 is therefore given by 

P(y a -)• v p ) = \{v p \v a (t))\ 2 

= E K.u^u;^-^. (7) 
k,j=i 

Here, we can approximate 



1 1 /\ Ttl^ 

E k - E 3 = ^ + ml- jj? + m) « (8) 
By setting t ps L, where L is the propagation length, we get 

3 A 2 

^) = E UakU pk U aj U; j e- i ^ L . (9) 

From the quantity appearing in the exponential one defines the oscillation phase 
2^ = ^IW~L = T~ w ^ n tne oscillation length L = For large values of the phase, 

the oscillations become more and more rapid and the observable effect is washed out. 
As a result, clear signals arise for small oscillation phases only. However, for too small 
oscillation phases, oscillations do not occur. Therefore, the oscillation length should be 
of the order of the propagation length. 

In several experimental regimes, the oscillation reduces to a two-flavour problem (for 
details, see [23]). In this case the basis change from flavour eigenbasis to mass eigenbasis 
can be characterized by a single angle 

u=( cose „ sia i). (10) 

\ — sin V cos v J 
The transition probability then simplifies to 

P(Ua _> Vp) = sin 2 (2 £) S in 2 {^^j , (11) 

where Am 2 is the difference of the squared masses of the two remaining mass eigen- 
states. From that we find the survival probability P{y a —> v a ) = 1 — P{y a — > vp)- 
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3. Neutrino oscillations with a minimal length 

Let us now implement the minimal length i. In conventional quantum mechanics, 
wave number and momentum coincide, k = p. From this, one finds that momentum 
eigenstates are plane waves in the position eigenbasis due to 



/ 



dk' e ik ' x 5(p(k) - p(k')) = J dp' e tp ' x 5(p - p') = e ipx , (12) 



where for simplicity we are considering the one-dimensional case. This situation changes 
in GUP models. As explained in detail in [24], the momentum remains unbounded in 
GUP theories. The wavelength of a particle, however, is restricted to be larger than the 
minimal length A > £, leading to a non-linear relation k(p) between the wave vector and 
the momentum. This is what also happens in the case of the coherent state approach to 
non- commutative geometry [10]. As mentioned above, all modifications due to the non- 
commutativity can be accounted for by modifying the momentum integration measure 
dp — > dpe~ £2p2 . This is equivalent to GUP theories for the choice f(p 2 ) = e £2p2 — 1 in 
equation (1). From this, we see that the Fourier transform of a momentum eigenstate 
is given by 



/ 



dk 



dp' 



ik(p')x ~ P' 
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(13) 



dke ikx 5{p{k)-p(k')) = J dp' 

= J dp' e -W e *V)* 5 b-J) = e *(p)* > (14) 
where 

fj* = e -<V ^ k{p) = ^ erf (£j9 ) + c (15 ) 

with the error function erf (a;). The integration constant c has to be zero, because k — 
when p — 0. This dispersion relation reduces to the classical case for p <^ i^ 1 and 
shows a saturation behaviour for large momenta with k(p) P =°° Generalising this 

to four-vectors k^ = (cu, ftj and p^ = (E,p) we have 



F = ^ erf (^ M )' ( 16 ) 
and in particular 

u{E) = ^erf(^). (17) 

Indeed, in [24] it is shown that the mass condition E 2 = p 2 + m 2 remains valid if u(E) 
has the same functional form as k(p). Starting from (17), we can follow the standard 
derivation of the oscillation probability (Eq. (7)- (9)). The crucial step is given explicitly 
by 
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^ (erfU\v*\) + — — e~^y 



Am 



2 



b ±i e -m\ (18) 

modified transition probability 



2£ 

where we have used ~ ~ i? and Am\- = ml — mj. From this we find the 
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Pt(v a = £ ^akUpkU^e-*^^- 1 E ' )L . (19) 

fej=i 

The two-flavour case can be obtained starting from equation (10). The modified 
transition probability reads 

P{y a ^) = |(^|^(t))| 2 = sin 2 (2^) sin 2 (^^y^) , (20) 

which on using approximation (18) and t = L gives 



P{v a ->■ ^) = sin 2 (26») sin 2 ^ 



Am 2 L 



iE e—J . (2!) 

As an illustration, we plot standard and modified oscillations for the two-flavour case 
in figure (1). 



4. Analysis of the modified oscillation pattern 

Let us now analyse the modified oscillation probability Eq.(19). Comparing our result to 
the literature, we find that in contrast to [15], the suppression factor for the oscillation 
phase is not dependent on the propagation length but only on the energy. 
The relative difference between the classical oscillations and the minimal length 
oscillations is shown in figure (2). To get a very rough estimate of the value of the 
minimal length parameter £, we apply our modified oscillation probability Eq.(19) to 
available data for — y transitions from the MINOS experiment. As in [25], we 
employ the two-flavour case Eq.(21) for the analysis using the standard oscillation 
prediction from [25]. The results are shown in figure (2). It is obvious that the 
differences between classical and modified behaviour are smaller than the experimental 
uncertainties for all realistic values of £. From a x 2 test of the data using the modified 
oscillation pattern we only find a very weak bound of l~ x > 10 GeV (95% CL). However, 
increasing the neutrino flux by a factor of 100 might allow to push the limit into the 
100 GeV region. Whether the size of the effect gets modified by a more sophisticated 
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Figure 1. Comparison of classical and minimal length oscillations for the two- 



flavour case with Am 2 
I- 1 = 10 GeV. 



2.3 • 10" 3 eV , sin"(26») = 0.9, baseline L = 5000 km and 



analysis is a subject of future studies. Thus, neutrino oscillations can, in principle, 
provide novel information on NCG-induced minimal length pheonomena. 

As an alternative, we look for experimental conditions which could lead to the 
detection of this phenomenon in the near future. To this purpose we test the parameter 
space for baseline and energy range of neutrino experiments by considering a significant 
departure from classical oscillation probabilities, i.e. 

Ap= \P(a->P) -Pi{a^P)\ > 0.1. (22) 



Table 1. Oscillation parameters used in the analysis of the modified oscillation 

behaviour. Values taken from [26]. 

Parameter Value Parameter Value 

sin 2 (26»i 2 ) 0.861 Am\ 2 7.59-10" 5 eV 2 

sin 2 (26>i 3 ) 0.15 Am 2 13 2.43-10" 3 eV 2 

sin 2 (26> 23 ) 0.92 Am 2 2?> 2.43-10" 3 eV 2 

Sep 



Fig. (3) shows regions in which the above condition for the experimental parameters 
of propagation length L and neutrino energy E for a — y transition in the three- 
flavour model is fulfilled, with oscillation parameters as in table (1). Measuring neutrino 
oscillations in these parameter regions should lead to an observable increase or decrease 
of detected muon neutrinos in comparison to the standard oscillations. Note that 
within these regions the number of neutrino counts can be larger by a factor of 10 or 
more for small values of the classical oscillation probability. As we are only looking at 
the difference in oscillation probabilities, these results are only weakly dependent on the 
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Figure 2. Upper part: Expected neutrino flux for the minimal length model with 
different fundamental scales from l~ x — 1 GeV to l~ x — 1 TeV (from top to bottom) 
for a V[i — y Vy, transition with baseline L ~ 734 km. MINOS data taken from [25]. 
Lower part: Relative differences in the oscillation probabilities for different 
fundamental scales. 



values of the oscillation parameters as further calculations confirm. We see that already 
at distances within the solar system a significant effect could be measurable. 

Looking at Eq.(19) we see that for neutrinos with energies above the fundamental 
scale oscillations are strongly suppressed compared to the classical case. This is 
equivalent to saying that the oscillation phase in the minimal length scenario, 

± = *£eM-e&), <*) 

is exponentially damped and coherence is maintained for much larger distances than in 
the standard oscillation framework. This is a clear experimental signature which could 
be tested using astrophysical sources of high-energetic neutrinos. Neutrinos of sufficient 
energy are for example created in the atmosphere by cosmic rays for which, however, the 
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Ap > 0.001 
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Figure 3. Parameter space for three-flavour oscillations and £~ x = 1 TcV. In the 
indicated regions neutrino oscillations are coherent and show a significant departure 
from classical oscillations for different values of Ap. 

baseline is too short to provide a realistic setup. The suppression of oscillations should 
also leave an imprint in the flux of high-energy cosmogenic neutrinos. As cosmogenic 
neutrinos are created in flavour eigenstates u e , only, neutrinos with energies greater 
than the fundamental scale will not oscillate, leading to a suppression of r-neutrinos in 
the cosmogenic flux. This, however, is difficult to establish experimentally due to the 
very small flux of cosmogenic neutrinos. 

On the other hand, possible point sources of high-energy neutrinos such as 
gamma-ray bursts or active galactic nuclei could provide the possibility to test the 
modified oscillation behaviour even without exact knowledge of the baseline and flavor 
composition of the neutrinos. In the standard oscillation scenario, whatever the 
composition of the high-energy neutrino beam at the source is, the neutrinos will reach 
the earth in a fully mixed state (if the baseline is larger than the oscillation length 
which we assume is the case for the aforementioned point sources). In contrast, in our 
minimal length scenario the high-energy neutrinos will reach the earth in their original 
flavor composition. If neutrino telescopes such as IceCube or ANTARES indeed find a 
deviation from a perfect flavor mixing for high-energy neutrino spectra coming from a 
point source, this clearly hints at a modification of neutrino oscillations such as the one 
advertised in this paper and allows a test of the proposed scenario within near future. 

5. Summary 

To summarise, we have implemented an effective quantum gravity minimal length. 
Modifications of the neutrino oscillation probability were calculated and show a 
deviation from the classical behaviour at high energies. A rough estimation using 
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MINOS data limits the size of the minimal length to (r 1 > 10 GeV. A significant 
increase in statistics of earthbound experiments or a drastic change of the experimental 
baseline will allow to put stringent limits on (r l . Therefore, we have considered 
alternative neutrino sources, such as atmospheric, galactic, extragalactic and cosmogenic 
neutrinos. We showed that gamma ray bursts and active galactic nuclei could provide 
signals of modified oscillation probabilities due to the presence of the minimal length. 
Neutrinos with energies larger than the fundamental scale would propagate from the 
source to the Earth in their original flavour state, with a complete suppression of 
oscillations. This is in contrast to the standard oscillation scenario and provides a 
striking experimental signature that would allow to test our current understanding of 
quantum gravity. 

This work is supported by the Helmholtz International Center for FAIR within the 
framework of the LOEWE program (Landesoffensive zur Entwicklung Wissenschaftlich- 
Okonomischer Exzellenz) launched by the State of Hesse. 
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